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Abstract. The extraordinary properties of laser-induced transparency of a 
negative-index slab and parametric amplification for a backward-wave signal are 
investigated. The effects of the idler absorption and phase mismatch on the 
amplification of the signal are studied, and the feasibility of ensuring robust 
transparency for a broad range of control field intensities and slab thicknesses is 
shown. A particular option consisting of the independent engineering of a strong 
four-wave mixing response and the negative refractive index is proposed and its 
specific features are investigated. The feasibility of quantum control over the 
slab transparency in such a scheme is confirmed through numerical experiments. 
We thus show opportunities and conditions for the compensation of the strong 
absorption inherent to plasmonic negative-index metamaterials, and we further 
show achievable transparency through coherent energy transfer from the control 
optical field to the negative-index signal. 
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Keywords: Nonlinear optics, Negative-index metamaterials, Coherent optical effects, 
Backward waves, Parametric microamplifiers and oscillators 



Coherent Nonlinear Optics and Quantum Control in Negative-Index Metamaterials 2 
1. Introduction 

Negative-index (also known as negative phase velocity or left-handed) metamaterials 
(NIMs) form a novel class of electromagnetic media that promises revolutionary 
breakthroughs in photonics [T]. Significant progress has been achieved recently in 
the design of bulk, multilayered, negative-index, plasmonic structures [21 [3]. The 
majority of NIMs realized to date consist of metal-dielectric nanostructures that have 
highly controllable magnetic and dielectric responses. The problem, however, is that 
these structures have losses that are difficult to avoid, especially in the visible range 
of frequencies. Irrespective of their origin, losses constitute a major hurdle to the 
practical realization of the unique optical applications of these structures. Therefore, 
developing efhcient loss-compensating techniques is of paramount importance. So 
far, the most common approaches to compensating losses in NIMs are related to 
the investigation of the possibility to embed amplifying centers in the host matrix 
[1]. The amplification is supposed to be provided through a population inversion 
between the levels of the embedded centers. Herein, we investigate alternative options 
based on coherent, nonlinear optical (NLO) energy transfer from the control optical 
field(s) to the signal through optical parametric amplification (OPA). Nonlinear optics 
in NIMs remains so far a less-developed branch of optics. On a fundamental level, 
the NLO response of nanostructured metamaterials is not completely understood or 
characterized and cannot be predicted effectively to date. Nevertheless, it is well 
established that local-field enhanced nonlinearities can be attributed to plasmonic 
nanostructures, and some rough estimates of their magnitude can be obtained. The 
feasibility of crafting NIMs with strong NLO responses in the optical wavelength 
range has been experimentally demonstrated in [4]. Unlike natural, positive-index 
(PI) materials, the energy fiow and the phase velocity are counter-directed in NIMs, 
which determines their extraordinary linear and NLO propagation properties. Unusual 
properties of nonlinear propagation processes in NIMs, such as second harmonic 
generation, three-wave mixing (TWM) and four-wave mixing (FWM) OPA, which 
are in a drastic contrast with their counterparts in natural materials, were shown in 
[SlelElilElITOliniinilinilTl]. striking changes in the properties of nonlinear pulse 
propagation and temporal solitons |15| . spatial solitons in systems with bistability 
[161 ITTl [18] , gap solitons [19] , and optical bistability in layered structures including 
NIMs [20l were revealed. A review of some of the corresponding theoretical approaches 
is given in [21]. 

In the present paper, we investigate the effects of idler absorption and phase 
mismatch on parametric amplification of the backward waves and propose a novel 
scheme of compensating losses based on the results of this investigation. The paper 
is organized as follows. In Section [2l the similarities between the solutions to the 
nonlinear propagation equations for the electric and magnetic, quadratic and cubic 
nonlinearities are shown for the case of uniform control fields. Here, the local linear 
and nonlinear parameters are assumed independent of the intensities of the control 
fields, and the conclusions are applicable to both TWM and FWM. The induced 
transparency exhibits a resonance behavior as a function of the control field intensity 
and the NIM slab's thickness due to the backwardness of the light waves in NIMs. 
Usually, the resonances are narrow, especially if the parametric process is assisted by 
amplification of the idler due to Raman or population-inversion gain. The sample 
remains opaque anywhere beyond the resonance magnitudes of the control field and 
the resonance thickness of the sample. Counterintuitively, we show that transparency 
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becomes achievable within a broad range of these parameters if the absorption for 
the idler exceeds that for the signal. Phase matching of the contrapropagating waves 
presents a technical challenge [HI [521 [H]. We also show that, in the above indicated 
case, the transparency of a NIM slab also becomes much more robust against the 
phase mismatch. Based on these outcomes, we propose and investigate in Section 
[3] the option of independent engineering of the negative index and FWM nonlinear 
response, which is different from the one proposed earlier [11 1 112 1 [13]. We consider the 
doped metamaterial, where the signal appears in the vicinity of the transition between 
the excited energy levels, whereas the idler couples with the ground state of the 
embedded quasiresonant centers that provide a resonantly enhanced, FWM response. 
Here, all local optical characteristics including nonlinear susceptibility exhibit a strong 
dependence on the intensity for the driving fields and on the frequency resonance 
offsets of the coupled waves. Hence, the output signal can be tailored through the 
means of quantum control. It is shown that the indicated cardinal changes in the 
coupling scheme bring about major changes in the properties of the laser-induced 
transparency of the doped NIM slabs. The possibility of eliminating the negative 
role of the phase mismatch on the tailored transparency of the slab in this case is 
also shown. The results, supported by numerical simulations, prove the feasibility 
of tailored transparency, amplification and the creation of a microscopic, mirrorless, 
backward-wave, optical parametric oscillator that generates contradirected beams of 
entangled right- and left-handed photons. 

2. Backward waves, parametric interaction in a NIM and solutions to the 
nonlinear propagation equations 

2.1. Poynting and wave-vectors in a lossless NIM 

We consider a traveling electromagnetic wave, 

E(r, t) = (l/2)Eo exp[i(k • r - ujt)] + c.c, (1) 
H(r, t) = (l/2)Ho exp[i(k • r - ut)] + c.c. (2) 

From the equations 

VxE = — — ,B = AiH,VxH=- — ,D = £E (3) 
c ot c at 

one finds that 

k X E = -/xH,k X H = --eE, V^£; = -^i/. (4) 
c c 

Equations ([4]) show that the vector triplet E, H, k forms a right-handed system for 
an ordinary medium with > and /i; > 0. Simultaneously negative and /x,; result 
in a left-handed triplet and negative refractive index 

n = -^, = n^{ujlcf. (5) 

We assume here that all indices of e, ^ and n are real numbers. The direction of the 
wave- vector k with respect to the energy flow (Poynting vector) depends on the signs 
of e and ^: 

S(r,t) = f [E X H] = = ^E^. (6) 

At Ei < and /i^ < 0, S and k become contradirected, which is in contrast with the 
electrodynamics of ordinary media and opens opportunities for many revolutionary 
breakthroughs. 
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2. 2. Coupling geometry and coherent energy transfer from the ordinary control fields 
to the backward signal in a NIM 



Figure 1. Coupling geometry for tliree-wave mixing (a), four-wave mixing (b) 
of the ordinary and one backward electromagnetic waves and scheme of quantum 
controlled four-wave mixing in the embedded resonant nonlinear-optical centers 
(c). Here, Si, ki and uii are energy flow, wave-vector and frequency of the 
negative-index signal (ni < 1); 83^4, 4 and aJ3^4 - of the positive-index control 
fields; k2 and aJ2 stand for the positive-index idler. 



The basic idea of compensating losses by coherent energy transfer from the control 
field(s) to the signal through parametric interaction is illustrated in Fig. [TJ We 
assume the wave at uji with the wave- vector ki directed along the z-axis is a negative- 
index (NI) (ni < 0) signal. Therefore, it is a backward wave because its energy flow 
Si ~ (c/47r)[Ei X Hi] appears directed against the z-axis. In the TWM case, shown 
in Fig. [ija), the medium is illuminated by a higher-frequency, ordinary PI wave at 
0J3 traveling along the z-axis (713 > 0). In the FWM case, Fig. [Ijb), the slab is 
illuminated by two PI control (pump) waves at 0^3 and oj^. In both cases, all wave- 
vectors are co-directed along the the z-axis. Due to the parametric interaction, the 
control and signal fields generate a difference-frequency idler at 0^2 = W3 — wi (TWM) 
or at UJ2 ^ UJ4 + UJ3 — LOi (FWM) , which is also assumed to be a PI wave (712 > 0) . The 
idler contributes back into the wave at wi through the same type of the parametric 
interaction and thus enables OPA at u>i by converting the energy of the control fields 
into the signal. Thus, all of the coupled waves have their wave- vectors co-directed 
along z, whereas the energy flow of the signal wave. Si, is counter-directed to the 
energy flows of all the other waves, which arc codirected with their wave-vectors. 
Such coupling schemes are in contrast both with the conventional phase-matching 
scheme for OPA in ordinary materials, where all energy-flows and phase velocities are 
co-directed, as weU as with TWM backward- wave mirrorless OPO [H [H mUS], 
where both the energy flow and wave-vector of one of the waves are opposite to all 
others. 

2.3. Equations for coupled contrapropagating backward and ordinary waves 

First, we shall show that magnetic and electric TWM and FWM processes can be 
treated identically. We consider two alternative types of nonlinearities - electric. 




(a) 



(b) 



(c) 
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D = eE + 4ttP^^, B = /iH; and magnetic, B = + AttM^^, D = eE. Nonlinear 
polarization and magnetization are sought in the form 

P^^(r, t) = (l/2)Po^^(r) exp[*(k • r - ijt)] + c.c, (7) 

M^^(r,t) = (l/2)Mo^^(r)exp[i(k-r-cjt)]+ c.c. (8) 

Accounting for Eqs. (|3]), one can derive 

V X V X E = -{fi/c^)d^B/dt'^, 

-AE = ^(c^Vc')[eE + 47rP^^], (9) 

V X V X H = -{e/c^)d^B/dt^, 

-AH = e(wVc')[MH + 47rM^^]. (10) 

For the medium with the electric nonlinearity, the equation for the slowly varying 
amplitude Eq of the wave with the wave-vector along the z-axis takes the form: 

dEo/dz = i^{2nuj^/kc^)P^^ exp[i(k - k)z]. (11) 
For the magnetic nonlinearity, the equation is 

dHo/dz = i£{2Tru^/kc^)M^^ cxp[i(k - k)z]. (12) 
The equations are symmetric and can be converted from one to the other by replacing 

/Li < — > €. 

For the electric quadratic nonlinearity. 

Pi"" - X^eiE^Ei cxp{i[(fc3 - k2)z ~ C^it]}, (13) 
= X^e2E3El exp{i[(fc3 - fcl)2 - U2t]}, (14) 

where uj2 = 1^3 — wi and kj = \nj\Lijj/c > 0. Then the equations for the slowly-varying 
amplitudes of the signal and idler in the lossy medium can be given in the form 

dEi/dz = i<JeiE^ exp[iAfcz] + {ai/2)Ei, (15) 
dE2/dz = i<Je2El exp[iAfcz] - {a2/2)E2. (16) 
(2) 

Here, aej = {kj /ej)2TrXej E3, Ak = ^3 — fc2 — fci, and aj are the absorption indices. 
The amplitude of the control (pump) wave E3 is assumed constant. 
For the magnetic type of quadratic nonlinearity, 

M^"- = X^\H3H; exp{z[(fc3 - k2)z - uj^t]}, (17) 

M^"" = X^nlHsH*, exp{z[(fc3 - h)z - c^2t]}, (18) 

the equations for the slowly-varying amplitudes are: 

dHi/dz = icTraiH^ exp[iAA:z] + {ai/2)Hi, (19) 

dH2/dz = ia,n2Hl exp[iAA:z] - (a2/2)i?2- (20) 

Here, dmj = (%/Mj)27i'X„y-ff3, Hz = const, and the other notations remain the same. 

For the electric-type FWM, the equations for the slowly-varying amplitudes are 
similar: 

dEi/dz =. ijiE^ exp[iAfc2] + {ai/2)Ei, (21) 
dE2/dz = ij2El exp[iAfcz] - {a2/2)E2. (22) 

Here, jj = {kj /ej)2T:xf^ E3E4 and Ak = k^ + k^ — ki — k2- 

We note the following three fundamental differences in equations P^ - lf^^ as 
compared with their counterpart in ordinary, PI materials. First, the signs of cti and 
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7i are opposite to those of (T2 and 72 because ei < and fii < 0. Second, the opposite 
sign appears with ai because the energy flow Si is against the z-axis. Third, the 
boundary conditions for the signal are defined at the opposite side of the sample as 
compared to the idler because the energy flows Si and S2 are counter-directed. 

We introduce effective amplitudes, ae,m,j, and nonlinear coupling parameters, 
ge,m,j, which for the electric and magnetic types of quadratic nonlinearity are defined 
as 

«ej = yJh/kj\E, , g,j = ^\k^k^/7^\2nxf^ E3, (23) 

amj = ^J\fij/kj\Hj,g„,j = V|/cifc2/Mi/^2|27rx|^]-i?3, (24) 
and for FWM as 

= \/\^3/kj\Ej , 5, = V\hh/^\2TTxf^ EsEi. (25) 

The quantities |ajp are proportional to the photon numbers in the energy fluxes. 
Equations for amplitudes aj are identical for all of the types of nonlinearities studied 
here: 

dai/dz= — giOj exp(iAA:z) + (Q!i/2)ai, (26) 
da2/dz = 52^1 exp(iAfcz) — (Q;2/2)a2. (27) 



2.4- Manley-Rowe relations and solutions to the equations for coupled 
counter-propagating waves 

At ai.2 = 0, gi = g2, e.g., for off-resonant coupling, one finds with the aid of equations 
© and (EH),' (EZl): 



d_ 

dz 



Slz 



S2Z 

hLU2 



0, 



a2 



0. 



(28) 



These equations represent the Manley-Rowe relations EH] , which describe the creation 
of pairs of entangled counter-propagating photons hioi and huj2. The second equation 
predicts that the sum of the terms proportional to the squared amplitudes of the signal 
and idler remains constant through the sample, which is due to the opposite signs of 
Slz and S2Z and is in contrast with the requirement that the difference of such terms 
is constant in the analogous case in ordinary nonlinear-optical materials. 

Taking into account the boundary conditions ai{z = L) = oil, and 02(2 = 0) = 
020 {L is the slab thickness), the solutions to equations (p6)) . ((27)) can be written as 

Afc 

ai{z) — Ai exp[(/3i + i-—)z] -f 



a^{z) 



A2 exp[(/32 + i 
KiAi exp[(/3i - 
K2A2 cxp[(/32 - 



2 

.Afc 

i—)z] 



2 
Afc 



where 



Pi,2 = (ai - a2)/4 ± iR, ki^2 = [±R + is]/g, 



R = Vr^, = 9*291, s = (ai + a2)/4 - i Afc/2, 



(29) 

(30) 

(31) 
(32) 
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Ak 

= {aiLK2 - a2oCxp[(/32 +i—)L]}/D, 

Ak 

M = -{fliL^i - cxp[(/3i + i — )L]}/D, 

Ak Ak 
D = K2 exp[(/3i +i — )L\ - ki exp[(/32 + i— 



(33) 
(34) 
(35) 



At Afc = and Im^ ~ 0, (ai + a2)L <JC tt (off- resonance) , equations ([^ and ([50)1 
reduce to 



al{z) 
a2{z) 



'■IL 



cos{gL) 



cos[gz) 



la 



IL 



sni{gz) + 



»a20 

cos{gL) ^ ' cos(gL) 

The output amplitudes are then given by 

^10 = [aiL/cos(,9L)] - ia2otan(5L), 
a2L = iai£tan(gL) + [a2o/cos(gL)]. 



sin[5(z - L)], 
cos[g(z- L)]. 



(36) 
(37) 

(38) 
(39) 



At 020 = 0, the equations for the energy distribution for the backward wave, 
Ti{z) = \ai{z 
take the form 



Ti{z) = |ai(z)/aii|^, and for the PI idler, ri2{z) — \a2{z) / ali^\'^ , across the slab 



Ti(z) = |[K2exp(/3iz) - Kiexp(/32z)]/-D| , 

|2 



(40) 

Tj2{z) = |[exp(/3iz)-exp(/32z)]/^|'. (41) 

Then the transmission factor for the backward- wave signal at z = 0, Tio, and the 
output idler aX z = L, r]2L, are given by 



Tw = 



V2L = 



ai(0) 



aiL 
02 (i) 



exp{-[(ai/2)-.s] L} 



aiL 



cos RL+ (s/R) sin RL 
{g/R) sin RL 



(42) 

COSRL+ (s/R) sin RL ' ^^^^ 

At aiL = 0, 02(2 = 0) = 020, the slab serves as an NLO mirror with a reflectivity 
(output conversion efBciency) at uii given by an equation identical to Ea. (|43)) : 



mo 



ai{0) 


2 


{g /R) sin RL ^ 


n* 
"20 




cos RL + (s /R) sin RL 



(44) 



3. Laser-induced transparency, amplification and generation of the 
backward wave 



The fundamental difference between the spatial distribution of the signal in ordinary 
and NI slabs is explicitly seen at aj = Ak = 0. There, equation (j42|) reduces to 

Tio = l/[cos(.gL)]2. (45) 

Equations ([5^ - (|i5)) show that the output signal and idler experience a sequence of 
geometrical resonances at gL (2j -I- l)7r/2, (j=0, 1, 2, ...), as functions of the 
slab thickness L and of the intensity of the control field (factor g). Such behavior 
is in drastic contrast with that in an ordinary medium, where the signal would grow 
exponentially as Ti cx exp(2gL). The resonances indicate that strong absorption of 
the left-handed wave and of the idler can be turned into transparency, amplification 
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and even into cavity-free self-oseillation when the denominator tends to zero. The 
conversion factors 7710 and 7720 experience a similar resonance increase. Self-oscillations 
would provide for the generation of entangled counter-propagating left-handed, fiwi, 
and right-handed, hLL}2, photons without a cavity. A similar behavior is characteristic 
for distributed-feedback lasers and is equivalent to a great extension of the NLO 
coupling length. It is known that even weak amplification per unit length may lead 
to lasing provided that the corresponding frequency coincides with high-quality cavity 
or feedback resonances. Below, we present investigation of the properties of such 
resonances. 

3.1. Effect of the idler absorption and phase mismatch on the laser-induced 
transparency resonances 

In order to demonstrate the major effects of the idler absorption and phase mismatch 
on the laser-induced transparency resonances, we consider the model where the 
dependence of the local optical and NLO parameters on the intensity of the control 
field can be neglected and the parameter g is real. Such a model is relevant to, e.g., 
off-resonant quadratic and cubic nonlinearities attributed to the structural elements of 
metal-dielectric nanocompositcs [3] . The results will be used in the next subsection for 
optimization of the transparency achievable through embedded resonant FWM centers 
with power-dependent optical parameters. A crucial role of the outlined geometrical 
resonances is illustrated in Figs. [2] and [3l Besides the factor g, the local NLO energy 
conversion rate for each of the waves is proportional to the amplitude of another 
coupled wave and depends on the phase mismatch Afc. Hence, the fact that the 
waves decay in opposite directions causes a specific, strong dependence of the entire 
propagation process and, consequently, of the transmission properties of the slab on 
the ratio of the decay rates [13] . A typical NIM slab absorbs about 90% of light at 
the frequencies which are in the NI frequency range. Such absorption corresponds to 
aiL « 2.3. Since the idler grows toward the back facet of the slab and the signal 
experiences absorption in the opposite direction, the maximum of the signal for the 
given parameters is located closer to the back facet of the slab. A change in the slab 
thickness or in the intensity of the control fields leads to significant changes in the 
distributions of the signal and idler along the slab and in their output values (Fig. [2]). 
As outlined above, the transparency exhibits an extraordinary resonance behavior as a 
function of the intensity of the control field and the NIM slab thickness, which occurs 
due to the backwardness of the light waves in NIMs and is illustrated in Fig. [3l 
Basically, such resonances are narrow, like those depicted in the plot corresponding 
to = 1 in Fig. [3{a), and the sample remains opaque anywhere beyond the 
resonance field and the parameters of the sample. If nonlinear susceptibility varies 
within the negative-index frequency domain, this translates into relatively narrow- 
band filtering. The slab would become transparent within the broad range of the slab 
thickness and the control field intensity if the transmission in all of the minimums is 
about or more than 1. Figure [3] shows the feasibility of achieving robust transparency 
and amplification in a NIM slab at the signal frequency through a wide range of the 
control field intensities and slab thicknesses by the appropriate adjustment of the 
absorption indices a2 > ai- Figure [3ja) depicts transmission properties of the NIM 
slab at aiL = 2.3 and different magnitudes of the absorption index > that 
are less than, equal to and greater than aiL. The figure shows dramatic changes 
in the transmission properties with changes in the ratio of the absorption indices 
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(a) (b) 

Figure 2. Tailored distribution of the signal, Ti{z) = \ai{z)/ai(L)\'^ , and of the 
idler, r]2(z) = \a2{z)/ai{L)\'^ , along the slab at aiL = 2.3, a2L = 3 in the first 
transmission minimum, (a): Afc = 0, gL = 3, 717. (b): Ak = I.Stt, gL = 4, 505. 

for the signal, ai, and the idler, a2- It is seen that the transmissions does not 
drop below 1 at a2 > cti- Hence, larger absorption for the idler is advantageous 
for robust transmission of the signal, which is counterintuitive. The increase of the 
idler's absorption is followed by the relatively small shift of the resonances to larger 
magnitudes of gL. Oscillation amplitudes grow sharply near the resonances, which 
indicates cavity-less generation. The distribution of the signal and the idler inside the 
slab would also dramatically change, as compared to that for the first transmission 
minimum depicted in Fig. [2l Unless optimized, the signal maximum inside the slab 
may appear much greater than its output value at ^ = 0. Phase-matching of the 
positive- and negative-index waves also presents a technical challenge. Figures[3jb)-(d) 
show the possibility to significantly diminish the negative role of the phase mismatch 
on the tailored transparency of the slab at the expense of a modest increase of the 
amplitude of the control field. At that, the spatial distribution of the signal and 
the idler may experience a dramatic change [Fig. [2] (b)]. Such dependencies are in 
strong contrast with their counterparts in PI materials and arc determined by the 
backwardness of the coupled waves that is inherent to NIMs. 

Only rough estimations can be made regarding x^^^ attributed to metal-dielectric 
nanostructures. Assuming x*^^^ ~ 10^^ ESU (~ lO'^pm/V), which is on the order of 
that for CdGeAs2 crystals, and a control field of / ~ 100 kW focused on a spot of 
D ~ bOfim in diameter, one can estimate that the typical threshold value of gL ^ 1 
can be achieved for a slab thickness in the microscopic range of L ~ l^m, which is 
comparable with that of the multilayer NIM samples fabricated to date [31 [5] . 

3.2. Independent engineering of NLO response, quasi-resonant four-wave mixing of 
ordinary and backward waves and quantum control of transparency of NIM slab 

Herein, we explore the feasibility of independently engineering the NI and the 
resonantly enhanced x'^'- response of a composite metamaterial with embedded NLO 
centers. We investigate resonant and quasi-resonant FWM and the accompanying 
processes that allow coherent energy transfer from the control fields to the counter- 
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(c) (d) 

Figure 3. Transmission of the negative-index slab, Ti{z), vs. parametric gain, 
gL, at aiL = 2.3 and different values of a2L and AkL. [(a), (c)] Ak = 0; [(d)] 
Ak = it; [(c), (d)] a2L = 4. 

propagating NI signal and PI idler. We show that this opens up opportunities 
for compensating optical losses in NIMs and for the creation of unique NIM-based 
photonic microdevices, which are of critical importance for the further development 
of nanophotonics in NIMs (for a review, see, for example, [T]). Among the specific 
features considered in this study are the above-outlined backwardness of the coupled 
waves and the possibility of quantum control over the nonlinear propagation process 
and its outcomes. The basic scheme of resonant four-wave mixing of the backward 
wave in a NIM is as follows. A slab of NIM is doped by four-level nonlinear centers 
[Fig. [IJc)] so that the signal frequency, uJi, falls in the NI domain, whereas all other 
frequencies, ws, lo^ and 0^2 , are in the the PI domain. Below, we show the possibility 
to produce transparency and even amplification above the oscillation threshold at 
uji controlled by two lasers at uj^ and UJ4 [Fig. [Hb)]. These fields generate an 
idler at 102 — uj^ + lo^ — uji^ which then contributes back in optical parametric 
amplification at wi. Unlike the scheme investigated in [TTl [121 US], here the idler 
corresponds to a higher-frequency transition from the ground state, and the signal 
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corresponds to a lowcr-frcqucncy transition between the excited states. No incoherent 
amplification is possible here for the idler, and the dependence of the idler and the 
signal absorption indices on the control fields cardinally changes. First, the scheme 
with relatively fast quantum coherence relaxation rates and the case where only a two- 
photon, Raman-like resonance for the signal holds is considered; all other one-photon 
frequency offsets are on the order of several tens of the optical transition widths. Then 
the scheme with the same quantum coherence relaxation rates but with higher partial 
spontaneous transition rates is considered, in which case population inversion at the 
coupled optical transitions is impossible. Finally we consider the scheme with longer 
quantum coherence lifetimes, which still does not allow population inversion at the 
optical transitions nor Raman-like amplification. The fact that all involved optical 
transitions are absorptive determines essentially different features of the overall loss- 
compensation technique in such composites in each proposed scheme. In all of the 
schemes outlined above, the linear and nonlinear local parameters can be tailored 
through quantum control by varying the intensities and frequency-resonance offsets 
for combinations of the two control driving fields. 

The following model, which is characteristic of ions and some molecules embedded 
in a solid host, has been adopted: energy level relaxation rates r„ = 20, Tg = Vrn = 
120; partial transition probabilities 7g„ = 50, 7mn = 70, (all in 10^ s~^); homogeneous 
transition half-widths F/g = 1, F/„i = 1.9, Tng = 1-5, F„,„ = 1.8 (all in 10^^ s~^); 
Tgm = 5, Tin = 1 (aU in 10^° s'^); Ai = 756 nm and Aa = 480 nm. The density- 
matrix method |28] is used for calculating the intensity-dependent local parameters 
while accounting for the quantum nonlinear interference effects. This allows us to 
investigate the changes in absorption, amplification, and refractive indices as well as 
in the magnitudes and signs of NLO susceptibilities caused by the control fields. These 
changes depend on the population redistribution over the coupled levels, which in turn 
strongly depends on the ratio of the partial transition probabilities. Electrical linear 
and nonlinear polarizations, Eq. ([7|), are calculated as 

Pi(z,<) = {l/2){P„\cxp{ikiz) 
+ exp[i(fc3 + ki — k2)z]} cxp(~iuJit) + c.c. 

= N{pngdgn + Pgndng); (46) 

+ Po2^ exp[i(fc3 + k^ — ^i)^]} exp(— iwai) 4- c.c. 

= N{pmldlm + Plmdml)- (47) 

Here, pij are the density matrix elements, and dij are the transition dipole elements. 

(3) 

Effective linear, Xi.2; and NLO, X121 susceptibilities dependent on the intensities of 
the driving control fields E3 and E2 are defined as 

Po\^XiEi, Po'i'^ = xf^ (48) 

Poi=X2i?2, Po'^'' = i?3i?4i?i*. (49) 
The linear susceptibilities determine the intensity-dependent contributions to 
absorption and to the refractive indices of the composite attributed to the embedded 
centers, while the NLO susceptibilities determine the FWM. Here, uji + UJ2 = + ^^4, 
and kj = \nj\iOjlc > 0. A significant difference between the resonant and off-resonant 
NLO processes is that all local optical parameters become intensity-dependent, and 
hence their spectral properties may experience a radical change near resonance. In 
particular, the NLO susceptibilities and, therefore, the parameters 71 and 72 become 
complex and differ from each other in the vicinity of the resonances. Hence, the factor 
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may become negative or complex. This indicates an additional phase shift between 
the NLO polarization and the generated wave that causes further radical changes in 
the nonlinear propagation features, which can be tailored. 

Figure [3] depicts such modifications at the given resonance offsets and intensities 
of the control fields. Here, Qi = uji ~ other resonance detunings Qj are defined 
in a similar way. Coupling Rabi frequencies are introduced as G3 = E3dig/2h and 
G4 = Eidnm/2h. The quantity 0:20 denotes the fully resonant value of absorption 
introduced by the embedded centers at 012 = uJmi with all driving fields turned off. 
Figure [D^a) displays the modified absorption/gain indices. The nonlinear spectral 
structures are caused by the modulation of the probability amplitudes, which exhibits 
itself as an effective splitting of the energy levels coupled with the driving fields. 
Figure Hl^b) shows the contribution to the phase mismatch associated with one such 
spectral structure. Figure Wic) and (d) indicate that the real and imaginary parts 
of the NLO susceptibilities become commensurate for the given susceptibility, but 
may exceed their counterparts for the idler by several times. This occurs due to 
the fact that different population differences contribute in different ways to the NLO 
susceptibilities [28], and driving fields cause significant redistributions of the level 
populations (Fig. [5]). At the given partial probabilities of spontaneous transition 
between the levels, population inversions at the signal transition become possible 
[Figl5]^e),(g)-(i)]. However, for the given frequency offsets of the control fields, 
corresponding amplification contributes negligibly to the energy conversion [Fig|4l^a)] . 
Alternatively, two-photon, Raman-like amplification at f2i ~ 20.05rg„ shown in 
FigHI^a) supports coherent, parametric energy-conversion from the control fields to 
the signal. Figures [6lja)-(d) display the spectral properties of the output signal 
at z = for one of the resonances in the vicinity of the signal frequency offset 
u>i — LUgn w 20.05rg„ at different optical densities of the slab at comi attributed to 
the impurity centers. We assume that the absorption of the host material in the slab 
at cui is fixed at 90% and it is equal to 88% at uj2- The density of the embedded 
centers and the slab thickness, and hence, the additional resonant optical thickness 
of the slab contributed by these impurities, may vary as shown in the panels. Actual 
quasi-resonant absorption/gain indices depend on the intensities and frequency offsets 
of the control fields, as shown in FiglJlJa). Besides the features imposed by the counter- 
propagation of the coupled waves, the output magnitudes of the signal at z = and 
the idler at z = L and their distributions inside the slab are determined by the 
interplay of several contributing linear and nonlinear processes. They include the 
phase mismatch, absorption of the signal and the idler, and the parametric gain g, 
which are all controlled by the driving fields E3 and E4. The dependence of the 
overall optimized output signal on the density of the impurities and on the slab 
thickness (on the resonant optical thickness of the slab) is depicted in Fig. Eljc). 
Such a behavior is determined by the radically different distributions of the idler, 
which propagates from left to right, and the signal, which propagates from right to 
left, [Fig. [6{f)]. Figures [6l^b)-(f) indicate the possibility of mirrorless self-oscillation. 
Figure [7] shows the role of partial spontaneous transitions between the energy levels. 
Here, 7,„„ = 9 x 10^ sec~^, which makes both population inversion and two-photon 
gain impossible [Fig. [Tlja)]. At the indicated Rabi frequencies and frequency offsets 
for the driving control fields, the energy-level populations are: n w 0.4, w 0.2009, 
Tn ~ 0.203f , w 0.2. The magnitude of the four- wave mixing coupling parameters 
appear comparable with those depicted in Fig. [4Uc)-(h). However, the absence of one- 
and two-photon amplification that would support energy-conversion processes, like in 



Coherent Nonlinear Optics and Quantum Control in Negative-Index MetamaterialslS 

[T^ and in Fig. |l{a), dramatically decreases the achievable amplification and increases 
the required optical thickness of the slab [Fig. [Itg),(h)]. Figure [5] shows that, even 
in such cases, the optimized magnitude of the required control field intensities and 
the slab optical density can be substantially reduced for centers with lower coherence 
relaxation rates and quasi-resonant coupling. Here, quantum nonlinear interference 
effects play an important role |28j . At the indicated Rabi frequencies and frequency 
offsets for the driving control fields shown in Fig. [8l the energy-level populations are: 
n w 0.504, rg w 0.165, r„ w 0.167, r„i « 0.164. Like in the previous examples, the 
losses in the host NIM material are taken to be fixed and equal to aNiMiL = 2.3 for 
the signal and aNiM2L = 2.1 for the idler. 

The dependencies presented in Figs. |4] - [8] correspond to the vicinity of 
the first geometrical resonance, which appears at the lowest magnitude of gL as 
shown in Fig. [3l For the above-indicated characteristics of the involved optical 
transitions, the magnitude G=150 GHz corresponds to control field intensities on 
the order of of / 10 kW/(0.1mm)^. Assuming a resonance absorption cross- 
section (740 ^ 10~^^ cm^, which is typical for transitions with oscillator strength 
of about one, and a concentration of embedded centers ~ 10^^ cm~^, we 
obtain a2o ^ 10^ cm~^ and the required slab thickness in the microscopic range 
L (1 — 100) fim. The contribution to the index of refraction by the impurities 
is estimated as An < 0.5(A/47r)a4o ~ 10^'^, which essentially does not change the 
negative refractive index. 

4. Conclusions 

The possibility to produce and to tailor a laser-induced optical transparency in a 
negative-index metamaterial slab through nonlinear-optical coherent energy transfer 
between ordinary control wave(s) and a negative-index backward signal is shown and 
proven by numerical simulations. Two possible types of nonlinear-optical couplings are 
discussed. One is off-resonant three or four-wave mixing by making use of the nonlinear 
susceptibilities that are assumed attributed to the metamaterial nanostructures and 
are independent of the intensities and frequencies of the coupled optical fields. The 
other option is the independent engineering of a resonantly enhanced, four-wave mixing 
nonlinearity associated with nonlinear-optical centers embedded in a negative-index 
host matrix. In the latter case, the proposed coupling scheme suggests that the 
frequency of the negative-index signal should fall in the vicinity of the transition 
between the excited levels of the centers, while the idler frequency appears coupled 
with the absorptive transition from the ground state. The scheme under investigation 
is different from the earlier proposed schemes and exhibits essentially different features. 
The extraordinary properties of the nonlinear-optical propagation processes in both 
outlined types of metamaterials are investigated. These properties are in drastic 
contrast with their counterparts in ordinary, positive index materials. The focus 
of this work is on the possibility of compensating for the strong losses inherent to 
metal-dielectric negative-index metamaterials and on producing laser-induced optical 
transparency and gain for the negative-index signal. In the case of a frequency- and 
intensity-independent nonlinear-optical response of the composite, the feasibility of 
producing transparency and amplification through the entire negative- index frequency 
domain above a certain control laser field intensity is investigated. This is shown to 
be possible by adjusting the absorption index for the idler to be greater than that for 
the negative-index signal. Specific features of the quantum control attributed to the 
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second scheme are investigated that aUow for transformable optics through frequency- 
tunable narrow-band transparency, quantum switching, filtering and amplification of 
light. The possibility of compensating the strong losses inherent to NIMs and realizing 
a miniature, mirrorless optical parametric generator of entangled contra-propagating 
backward and ordinary waves is also shown and supported by numerical simulations. 
The extraordinary features predicted in this work stem from the backwardness of 
electromagnetic waves, which is a feature inherent to this type of metamaterial. 
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Figure 4. Nonlinear spectral structures in local optical quantities produced by 
the control fields. y\ = {u)\ — uign) gn, ^^2 = aJ3 + 0^4 — . (a): absorption/gain 
indices for the signal and tfie idler; (b): phase mismatch; (c)-(f): four-wave mixing 
coupling parameters. Coupling Rabi frequencies and resonance frequency offsets 
for the control fields are: G3 = 254.64 GHz, ^3 = SOTgi, G4 = 108.48 GHz, 

S^4 = 20r77T,7i • 
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Figure 5. DifTcroncc of the energy- level populations and their dependence on 
the Rabi frequency of the control fields G3 and G4 (given in GHz). £^3 = SOFg;, 
= 20rm„. (h): G3=254.64 GHz, (i): G4=108.48 GHz. 
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Figure 6. Dependence of the transmission of the slab on the resonance frequency 
offset y\ = (oJi — ujgn)/Tgn for different optical densities of the slab, (a)-(d), on 
the resonant optical density of the slab, (e), and the distribution of the signal 
and the idler along the slab, (f). G3=254.64 GHz, Q3 = SOFg;, G4=108.48 GHz, 

1^4 = 2,0T^ mn ■ 
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Figure 7. Energy-conversion in the scheme with neither population inversion nor 
two-photon gain possible [fmn = 9 X 10^ sec~^); all other relaxation parameters 
are the same as in the previous case] . yi = {LUi—uign)/Tgn, ^^2 = t<J3-|-aJ4— tJi. (a): 
absorption indices for the signal and the idler; (b): phase mismatch; (c)-(f): four- 
wave mixing coupling parameters; (g) and (h): transmission factor, the dashed line 
shows transmission at g = 0. Coupling Rabi frequencies and resonance frequency 
offsets for the control fields are: G3=37.5 GHz, G4=202 GHz, = -Q^ = SOTgi. 
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Figure 8. Quasi-resonant coupling at lower quantum coherence relaxation 
rates and at neither population inversion nor two-photon gain possible [rg;=1.8, 

r,„„=i.9, rg„=i, r„,=i.5, r^g = 5 x 10-2, = 5 x 10-3 (in 10" 

sec^^); all other relaxation parameters are the same as in the previous case]. 

)/Tgn, u}2 = UJ3+LO4-UJ1. (a): absorption indices for the signal and 
the idler; (b): phase mismatch; (c)-(e): four- wave mixing coupling parameters; 
(f): energy level populations; (g) and (h): transmission factor, dash line shows 
transmission at g = 0. Coupling Rabi frequencies and resonance frequency offsets 
for the control fields arc: G3=0.86 GHz, ^3 = STgi; G4=8.99 GHz, Q4 = -2r™,„. 
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